A Unified Description of Screened Modified Gravity 
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We consider modified gravity models driven by a scalar field whose effects are screened in high 
density regions due to the presence of non-linearities in its interaction potential and/or its coupling 
to matter. Our approach covers chameleon, f(R) gravity, dilaton and symmetron models and allows 
a unified description of all these theories. We find that the dynamics of modified gravity are entirely 
captured by the time variation of the scalar field mass and its coupling to matter evaluated at the 
cosmological minimum of its effective potential, where the scalar field sits since an epoch prior to 
Big Bang Nucleosynthesis. This new parameterisation of modified gravity allows one to reconstruct 
the potential and coupling to matter and therefore to analyse the full dynamics of the models, from 
the scale dependent growth of structures at the linear level to non-linear effects requiring iV-body 
simulations. This procedure is illustrated with explicit examples of reconstruction for chameleon, 
dilaton, f(R) and symmetron models. 
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I. INTRODUCTION 

The discovery of the acceleration of the expansion of 
the Universe [1] has led to a reappraisal of some of the 
tenets of modern cosmology. In particular, the possibility 
of modifying the laws of gravity on short or large scales 
is taken more and more seriously [2]. 

In view of Weinberg's theorem stating that any Lorentz 
invariant field theory involving spin- 2 fields must reduce 
to General Relativity (GR) at low energy [3], any attempt 
to modify GR must involve extra degree(s) of freedom. 
The majority of known models involve scalar fields and 
can be separated into two broad classes, the ones involv- 
ing non-linearities in the kinetic terms and others with 
non-linear interaction potentials. All these models have a 
coupling of the scalar field to matter and there could be 
an environmental dependence which would manifest itself 
in the screening behaviour of the scalar field in high den- 
sity regions [4, 5]. Examples of such models abound: the 
dilatonic models [6, 7] generalising the Damour-Polyakov 
mechanism [8] where the coupling to gravity turns off in 
dense environments, the chameleon models [9, 10, 12-14] 
where a thin shell shielding the scalar field in dense bod- 
ies is present, the symmetron models [15-21] where the 
scalar field has a symmetry breaking potential where the 
field is decoupled at high density. 

Some models are essentially spin-offs of the previous 
ones like the f(R) theories [22-32] (for recent reviews 
of the f(R) gravity see [33, 34]) which are only valid 
when they behave like chameleon theories with a thin 
shell mechanism in dense environments [32]. In all these 
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examples, the large scale properties on cosmological dis- 
tances are intimately linked to the small scale physics as 
probed in the solar system or laboratory tests of grav- 
ity. Stringent constraints on the possible modifications of 
gravity follow from the cosmology of these models too. In 
particular, they may lead to potentially lethal variations 
of particle masses or Newton's constant during Big bang 
Nucleosynthesis (BBN). This must be avoided at all cost 
as this may destroy the formation of elements, one of the 
big successes of the Big Bang model. Such a catastrophe 
can be avoided provided the scalar fields sit at the min- 
imum of the density dependent effective potential prior 
to BBN. If this is the case, then the minimum of these 
models is stable enough to prevent large excursions of 
the scalar field and therefore of scalar masses/Newton's 
constant when the electron decouples during BBN. One 
of the most important consequences of this fact, which is 
common to chameleons, dilatons and symmetrons is that 
the cosmological background with the scalar field at the 
density dependent minimum of the effective potential be- 
haves essentially like the A-Cold Dark Matter (ACDM) 
model and is therefore almost indistinguishable from a 
cosmology comprising matter, radiation and a pure cos- 
mological constant. This is a major drawback and would 
immediately render irrelevant the modified gravity/dark 
energy models with screening properties. 

Fortunately, this is far from being the case as first an- 
ticipated in [10, 11] where the equation governing the 
density contrast of CDM was first studied. Indeed, in- 
side the Compton wavelength of the scalar field, the den- 
sity contrast grows anomalously compared to its usual 
growth in the matter dominated era. If this discrepancy 
were large enough on astrophysical scales, this may be 
detectable by future galaxy surveys. It turns out that 
the perturbation equation at the linear level depends on 
the time evolution of the scalar field mass and the cou- 
pling strength to matter. With these two functions, all 
the time and space properties of the linear perturbations 



can be calculated. 

In fact, these two time-dependent functions capture a 
lot more about the modified gravity models with screen- 
ing properties: they allow one to reconstruct fully and 
uniquely the whole non-linear dynamics of the models 
[5, 35]. Hence given these two functions, not only can 
one compute linear perturbations, but one can study the 
gravitational properties of the models in the solar sys- 
tem and laboratory experiments. One can also analyse 
the cosmological behaviour of the models with TV-body 
simulations. This way of defining the models, a reversed 
engineering procedure from the mass and coupling func- 
tions to the non-linear dynamics, is a lot more versatile 
than the usual direct route where a model is defined by 
its Lagrangian comprising the kinetic terms and an in- 
teracting potential. Indeed, all the usual models such as 
chameleons, /(-R), dilaton and symmctrons can be ex- 
plicitly rediscovered by specifying the particular ways the 
mass and coupling functions behave in time. Moreover, 
one can design new families of models. At the linear level 
of cosmological perturbations, this approach is equivalent 
to a space and time dependent parameterisation [36-45] 
in terms of the two Newtonian potentials obtained in the 
Jordan frame: the modified Poisson equation and the con- 
stitutive relation linking the two Newtonian potentials 
are directly and uniquely determined by the mass and 
coupling functions in the Einstein frame. For instance, 
we shall see below that one recovers the phenomenolog- 
ical description of f(R) models which uses a space and 
time dependent parameterisation [40] as a simple appli- 
cation of our formalism. 

The paper is arranged as follows, in a first part we de- 
scribe the modified gravity models with scalar fields and 
their cosmological background and gravitational proper- 
ties. We then describe the tomography of the models, 
i.e. how to reconstruct their full dynamics using the time 
evolution of the mass and coupling functions. In section 
IV, we focus on f(R) models. In section V we analyse the 
growth of structure. In section VI, we consider the con- 
straints on these models resulting from the variation of 
the fundamental constants. We conclude in section VII. 

Throughout this paper the metric convention is chosen 
as (— , +, +, +); Greek indices (fi,i/,- ■■) run over 0, 1, 2, 3 
while Latin indices (i,j, k, ■ ■ • ) run over 1, 2, 3. We shall 
adopt the unit c = 1 and mpi denotes the Planck mass. 
Unless otherwise stated a subscript o will always mean 
the present-day value of a quantity. 



the interaction potential together with the coupling of 
the scalar field to matter can be reconstructed from the 
knowledge of the scalar field mass m(a) and the coupling 
strength /3(a) as functions of time when the field sits at 
the minimum of the density dependent effective poten- 
tial. 

This mechanism relies on the fact that the scalar field 
must track that minimum since before BBN in order to 
preserve the constancy of particle masses at this epoch. 
In this section, we recall the setting of scalar field models 
and analyse their background evolution. 



A. Modifying Gravity with a Scalar Field 

The action governing the dynamics of a scalar field 
in a scalar-tensor theory is of the general form 



S = 




^R-\(v^ 2 -vW 



g£ rn (ip%,g^ 



(1) 



where g is the determinant of the metric g^„, R is the 

Ricci scalar and ipm are various matter fields labelled by 
i. A key ingredient of the model is the conformal coupling 
of (j) with matter particles. More precisely, the excitations 
of each matter field ipm couple to a metric g^ v which is 
related to the Einstein-frame metric g M „ by the conformal 
rescaling 



9tiu = A 2 (cj))g 



liv 



(2) 



The metric g^ u is the Jordan frame metric. Wc will anal- 
yse these models in the Einstein frame and come back to 
the Jordan frame picture later. 

The fact that the scalar field couples to matter implies 
that the scalar field equation becomes density dependent. 
More precisely, the scalar field equation of motion (EOM) 
is modified due to the coupling of the scalar field </> to 
matter: 



□0= -/3T- 



dV 



(3) 



where T is the trace of the energy momentum tensor T^, 
□ = V^V^ and the coupling of (j> to matter is defined by 



II. MODIFIED GRAVITY 



dlnA 



mm- 



(4) 



In this paper we propose a parameterisation of a broad 
class of theories with a scalar degree of freedom, such 
as the chameleon, dilaton and symmetron theories, and 
f(R) gravity. The success of these theories relies on mech- 
anisms that suppress the fifth force in local, high matter- 
density, environments. Wc will find that the complete 
non-linear Lagrangian comprising the kinetic terms and 



This is equivalent to the usual scalar field EOM with the 
effective potential 



V eS {4>)=V{<f>)-[A{4>)-l]T. 



(5) 



The role of this effective potential V r (</>) is crucial in all 
the modified gravity models we will consider. In essence, 



the effective potential is required to possess a unique mat- 
ter dependent minimum in the presence of pressure-less 
matter where T = —p m . The resulting potential 



has a minimum 
the minimum 



V oS {<S>) = V{cj>) + [A(<p) - l]p m (6) 

(p m ). The mass of the scalar field at 
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d 2 V cS . 



(7) 



must be positive. In many cases (such as the generalised 
chameleon and dilaton models discussed below) V{<j>) is 
a decreasing function and /3(</>) is an increasing function 
as <j>, though this is not the case for the generalised sym- 
metron model 1 . This guarantees that the effective poten- 
tial always has a minimum. In a cosmological setting we 
will also impose that m 2 3> H 2 with H being the Hubble 
expansion rate. It can be shown easily that, depending on 
the shapes of V{4>) and f3(4>), the chameleon, f(R), dila- 
ton and symmetron models are all described in a such a 
way. 

When matter is described by a pressure-less fluid with 



T^ = p m u»u v 



(8) 



and u^ = dx^/dT where r is the proper time, the matter 
density p m is conserved 







(9) 



where 6 = V^w^ and the trajectories are determined by 
the modified geodesies 



6 d^6 

u^ + ^—u" = -/?—. 
mpi mpi 



(10) 



In the weak-field limit with 

ds 2 = -(1 + 2$Ar)dt 2 + (1 - 2$ N )dx i dx i , (11) 

and in the non-relativistic case, this reduces to the mod- 
ified geodesic equation for matter particles 



dt 2 



-d l (<S> N + lnA((J))) 



(12) 



This can be interpreted as the motion of a particle in the 
effective gravitational potential defined as 



* = <I>Ar + hi ,4(0), 



(13) 



and is clearly a manifestation of the dynamics of modified 
gravity. 



1 For the generalised symmetron models, the potential is not mono- 
tonic but has the shape of a Mexican hat. However, in the part 
of the potential which will be of interest here, it is monotonically 
decreasing. 



When a particle of mass M in a homogeneous back- 
ground matter density is the source of gravity, the scalar 
field satisfies 



mpi 



(14) 



in which 5^ (r) is the 3-dimensional Dirac S- function and 
m the scalar field mass in the background, implying that 



* = -(! 



2(3 2 e- mr ) ^-, 



(15) 



where Gm = 



)~ 1 m,p 2 is the Newton constant. When 
(3 ~ O(l) and rrT 1 ^> r, this implies a substantial de- 
viation from Newton's law. For bodies much bigger than 
a point particle following the modified geodesies, non- 
linear effects imply that the effective coupling felt by the 
body is much smaller than (3 or the mass becomes much 
larger than the inverse of the typical size of the body 
(?7i _1 <C r). This is what happens in the chameleon model 
and f{R) gravity (the latter) and the dilaton and sym- 
metron models (the former), and guarantees that solar 
system and laboratory tests of gravity are evaded. 



B. Screening of Modified Gravity 

In this section, we shall unify the description for the 
screening 2 mechanisms [4, 5] involved in the chameleon, 
f(R) gravity, dilaton and symmetron models. As we shall 
see, the screening of large and dense bodies can be ex- 
pressed with a single criterion generalising the thin-shell 
condition for the chameleon models. The constraints we 
find are typically stated in terms of the scalar field mass 
mo in the cosmological background today and the cur- 
rent Hubble scale Hq, making £ = Ho/rriQ a key quantity. 
Physically, £ represents the range of the scalar fifth-force 
to the Hubble radius and a particular value that will be 
recurrent is mo/Ho ~ 10 3 or £ ~ 10~ 3 . This value means 
that the scalar field leaves its mark up to scales of the 
order of mega parsec, which again signals the transition 
where the modifications of gravity can be seen on linear 
perturbations or not. 



1. Chameleons 

The chameleon models (at least in their original form 
[9], see [10, 12-14] for other proposals) are characterised 
by a runaway potential and a nearly constant coupling 
(3. Chameleons are screened deep inside a massive body, 



2 To be clear, the 'screening' of a body refers to the fact that the 
deviation from Newtonian gravity, i.e., the fifth force exerted 
by this body on a nearby test mass, is suppressed to evade local 
constraints - in analogy to the screening of the electric force from 
a charged particle. 



where the field settles at the minimum </> c of V B s(4>) and 
stays constant up until a radius R s close to the radius of 
the body, R. In this case, the field profile is given by 

4> = 4>c R<R S (16) 

The field varies sharply inside a thin shell according to 

Pi. 



j__d 

r 2 dr dr 

and decays outside 



0- 



mpi 



R s <r < R 



Attiti 



pi 



1- 



R 



R> 



M e -m^{r-R) 



(17) 



(18) 



where </>oo is the minimum of the effective potential out- 
side the body and m^Af are respectively the masses 
of the scalar field and the body. At short distance com- 
pared to the large range m^ 1 , the effective gravitational 
potential is 



* = /3^L 



G N M 



mpi 



l -% 



(19) 



Gravity is strongly modified by a factor (1 + 2/3 2 ) if there 
is no shell inside the body (i.e., R s = 0) and one retrieves 
GR when i?, is close to R where 



AR 



\4>oo - <t>c\ 

6/3mpi$ 



N 



(20) 



with AR = R — R s and $_/v = GnM/R is the Newtonian 
potential at the surface of the body. The mass is screened 
when 



|0oo-0c| <2/3mpi$ 



N, 



(21) 



which is also the criterion to have a thin shell. 

More precisely, this implies several very stringent ex- 
perimental constraints on the chameleon models. The 
first one comes from the Lunar Ranging experiment [49] 
which measures the acceleration difference between the 
Earth and the Moon in the gravitational field of the Sun 



V 



^V&earth ^moon ) 
^carth i Ctmoon 



< io 



-13 



For the chameleon model we have [9] 



i?«j8- 



(AR. 



implying that 



P 



AR 



R- 



\ Rq 

l < io- 7 



(22) 



(23) 



(24) 



The Cassini experiment [50] imposes that the modifica- 
tion of the unscreened Cassini satellite in the vicinity of 
the sun should be such that 



/3 : 



ARc 

R(7) 



< 10" 



(25) 



Another type of constraint comes from cavity experi- 
ments where two small test bodies interact in a vacuum 
cavity [51]. This implies that 



P 



AR C 

Rca 



< 10 



-3 



(26) 



Finally, a loose bound must be imposed to guarantee that 
galaxies are not far off from being Newtonian [52] 



AR 



R 



Sal < ]_ 



gal 



(27) 



otherwise the modifications of gravity would have been 
seen by now in observations of galaxy clusters. These 
constraints strongly restrict the parameter space of the 
chameleon models. 



2. Symmetrons 

Symmetrons [17-21] are models with a mexican hat 
potential, a local maximum at the origin and two global 
minima at ±</>* like for example 



V(<t>) = V 



v 2 4>l 



(28) 



In general the term (</>/0j.) 4 can be replaced by any even 
function which is bounded below, without changing the 
qualitative properties of the model. 
Meanwhile, the coupling behaves like 



A{4>) = 1 + -f 



M ,2 



(29) 



close to 4> = 0. 

Let us consider a spherically dense body that is embed- 
ded in a homogeneous background. Inside this body the 
matter density p m is constant and the scalar field profile 

is 



C^sinhm^, r<R, 



(30) 



where the scalar field mass is given by m 2 = A-2p m — fJ 2 
and — [x 1 is the negative curvature of the potential V(<f>) 
at the origin. The field outside the body, on scales shorter 
than the large range m^ 1 associated to the scalar field 
value 4>oa which minimises V c fi((f)) outside, is 



D 



r > R, 



(31) 



where 



C = 
D = 



m c cosh™ CJ R' 
sinh m c i? — m c R cosh m c R 
m c cosh m c i? 



(32) 



If the body is dense enough, we have m 2 ps A 2 p m and 
m c R 3> 1, implying that D ps —Rcp^. Identifying the 
coupling to matter /3oo = mpiA24> QO , we find that the 
modified Newtonian potential outside the body is 



vp 



G N M 

r 
G N M 



1 + 



A 2 



1 







°4: 



A 2 m|i$jv 



■0(5: 



(33) 



for r sufficiently large compared to R. For R <C r -C m x ' 
the fifth-force is screened provided 



For a generic symmetron potential we have 4 (p 2 x ~ £^<p 2 

where 0+ is the minimum of T4g(0) in the cosmological 
background with matter density p^. Using ■ £2£ - ~ 10~ 6 , 
this leads to 



6 Si 



io-°^ 



1 



2A 2 m2, 1 $ Q 



1(1-6 rr2 

i^i£o<10-5 (39) 



which is easily satisfied for mo/Ho ~ 10 3 . Finally, in 
cavity experiments, the field inside the cavity is almost 
identical to the field in the bore, i.e., ~ 0, implying no 
deviation from usual gravity in such experiments. 



2^2771^$^ > 1, 



which is equivalent to 

|0oo - < 



< 2m P i/3 00 $ 



N- 



(34) 



(35) 



where <p c — 0. Note that this is the same screening crite- 
rion as in the chameleon case. 

The screening in the symmetron model depends on A 2 , 
$at and the environment through the environmental field 
value 0OO . Two test masses which are not screened when 
put in vacuum will be screened by a factor ((p^/cp*) 2 if 
they are in a region of high matter density (which implies 

0oo < 0*)- 

The transition of the minimum of V ff(<p) from = 
to = (p^ in the cosmological background happens in the 
recent past of the Universe provided 



/'" 



A 2 p 



'mOt 



where p m o is the present matter density. For a polynomial 
potential V(0), the mass-squared m 2 at the minimum 0* 
is of order p 2 , implying that the mass of symmetrons in 
the present cosmological background satisfies 



A 2 m Pl H Q , 



(37) 



One may see effects of modified gravity on astrophysical 
scales when tuq/Hq < 10 3 which implies that A 2 m,p l < 
10 6 . 

Using the screening criterion we find that the Sun and 
the Milky Way with $ Q ~ 10~ 6 arc marginally screened 
whereas the Earth with $0 ~ 10 -9 and the Moon with 
$ moon ~ 10~ n arc not screened. However, for the solar 
system tests such as the Lunar Ranging experiment 3 and 
the Cassini satellite, what is more relevant is the value of 
the symmetron field ga i in the Milky Way, which deter- 
mines the strength /3(0 ga i) of the modification of gravity. 

This imposes 



A2P 2 



al 



$.: 



<10" 



(38) 



3. Dilaton 

Dilatonic theories [6, 7] are very similar to symmetrons 
in as much as they share the same type of coupling func- 
tion, 



A(0) = l + #( 



M 2 , 



(40) 



but they differ as the dilaton potential V{<j>) is a mono- 
tonically decreasing function of 0. All the dynamics can 
be analysed in the vicinity of 0* as the minimum of the 
effective potential is close to 0* for large enough A 2 . 
The density dependent minimum of V e s{(p) is given by 



0min(/Om) ~ 0* = 



A 2 p m ' 



(36) with the mass given by 



A 2 p„ 



(41) 



(42) 



where m* = 771(0*) and the potential is chosen to be a 
quintessence potential such that m 2 ~ Hq. 

Let us consider a spherically dense body. Inside the 
body we have 



and outside 



sinh m c r _ 

= C + C* — , r<R, 



D 

3 



(43) 



(44) 



for distances shorter than the range m^ 1 . When m c R 3> 
1, we find that 



D PS -R(4>oo -0 C ), 
and the effective Newtonian potential is 



* = — 



G N M 



1 



^2(0oo - 0c)(0oo ~ 0*) 



<l> 



N 



(45) 



□2 

+ 0(-s-), 



(46) 



The Nordtvcdt effect leads to a weak bound [17]. 



See Eq. (19) in [17] for a more accurate expression. 



for R <C r <C m^ . Outside the body we have 

A2mpi 

with /Soc = /3(4>oo) an d therefore 

mpi 



from which we deduce that 



and finally 



A 2 mpi 

B 2 
A 2 ml l <£> N 



Poc 

Pc 



Poc 

Pc 



(47) 



(48) 



(49) 



o(5)- 



(50) 

for i? -C r <C JTIqc , . The screening criterion is (almost) 
the same as in the symmetron case 



2A 2 TOpi$AT > I 1 - 



Pc 



or equivalcntly 

|0o 



C | < 2/3(</) 00 )mpi$7v, 



(51) 



(52) 



which is the same as in the chameleon and dilaton cases. 
The mass of the dilaton today in the cosmological back- 
ground is 



A2Pmo = 3yl2m P1 f2 m oi?o > 



(53) 



in which f2 m o is the present value of the fractional energy 
density of matter Q m , implying that ^mpj ~ 10 6 for 
models with mo/Ho ~ 10 3 . 

As in the symmetron case, this implies that both the 
sun and the Milky Way are marginally screened when 
surrounded by the cosmological vacuum. But given that 
what matters for the magnitude of modified gravity is the 
dilaton value t^oo = ga i in the Milky Way, the Cassini 
bound can be written as 



A, 



gal - 0c)(0gal - </>*) 



$ 



< 10 



-5 



N 



which leads to 



A 2 m pl $ Q Pgal 



< 10" 



(54) 



(55) 



Using ■ £22 - ~ 10 6 , we see that the Cassini bound is sat- 

Pgal 

isficd for dilatons. 



where <f> c is the value inside the body assumed to be at 
the minimum of the effective potential, (j)^ the minimum 
value outside the body and <J>at is Newton's potential at 
the surface of the body. This is a universal criterion which 
is independent of the details of the model. In fact, it de- 
pends only on the values of the scalar field which min- 
imises the effective potential V s(<f>) inside and outside 
the body. If this criterion is satisfied, then the value in- 
side the body does not deviate much from the minimum 
value there. 

Phenomenologically, we have just recalled that strin- 
gent local constraints on modified gravity can be ex- 
pressed in terms of the screening condition. In the fol- 
lowing we shall assume that the Milky Way satisfies the 
screening criterion. When this is the case, local tests of 
gravity in the solar system and in the laboratory can be 
easily analysed as </> ga i can be determined analytically. 
In the chameleon, dilaton and symmetron cases, this al- 
lows one to determine bounds on the ratio ttiq/Hq which 
essentially dictates if modified gravity has effects on as- 
trophysical scales. The screening condition for the Milky 
Way may be relaxed slightly for some model parameters 
because it is itself in a cluster with higher density than 
the background. In this case, full numerical simulations 
are required to determine ga i and see if local tests of 
gravity are satisfied. This may enlarge the allowed pa- 
rameter space of the models slightly and lead to inter- 
esting effects. Numerical simulations are left for future 
work. 

One of the advantages of the screening condition is that 
it only depends on the minimum values of the scalar field 
in different matter densities. In the following section, we 
will find an explicit formula for <p c — 0oo which depends 
only on the time variation of the mass 771(a) and coupling 
/3(a) in a cosmological background. This may seem sur- 
prising as the behaviour of the scalar field may appear 
to be loosely connected to the scalar field dynamics in 
a static environment. In fact, the relation between both 
regimes of modified gravity, cosmological and static, fol- 
lows from the fact that the scalar field sits at the mini- 
mum of its effective potential V c s(4>) since before BBN. 
As it evolves from BBN through the dark ages and then 
the present epoch, the cosmological values of the scalar 
field experience all the possible minima of V e s ((/>). Hence 
realising a tomography of the cosmological behaviour of 
the scalar field, i.e., just knowing its mass and coupling 
to matter as a function of time since before BBN, will 
allow us to analyse the gravitational properties of the 
models. 



4- The Screening Criterion 

We have seen that all the models of the chameleon, 
dilaton and symmetron types lead to a screening mecha- 
nism provided that 



|0oo - 0c I < 2/3(0 oo )m P i<I>7V, 



(56) 



5. The Reason for a Universal Screening Condition 

As we have seen in the examples above, we get the 
same screening condition for all known models. Below 
we argue why this is the case for a whole range of models 
satisfying only some simple assumptions. 

We start with the most general model for the behaviour 



of the scalar field in matter 



V ei 



V, 



P(4>)Pn 

mpi 



(57) 



and we will analyse the standard setup - a spherical body 
of density p c and radius R embedded in a background of 
density p^ - under the following assumptions: 

1 . The effective potential has a matter dependent min- 
imum 4>{p). 

2. For any (physical) solution to the field equation, 
the mass of the field at r — 0, ms = m(cf>s, p c ), is 
a positive monotonically increasing function of the 
density p c and satisfies 5 lim Pc _ >00 m(<fis(pc),Pc) = 
oo. 

3. Outside the body, where p^ <C p c , and within 
the Compton wavelength of the field ra^ 1 the solu- 
tion to the field equation is well approximated by 
4> = 0oo + — • This means that a first order Taylor 
expansion around (f)^ holds outside the body. 

Now we can look at the solutions to the field equation 
under the previous assumptions. The field starts out at 
some field-value </> = 4>s inside the body, and close to 
r = the solution can therefore be written 

/sinh(msr) 



<f> = <j> s + B[ ^-^ - f (58) 

V m sr J 

for some constant B. We can for our purposes, without 
loss of generality, assume that B > 0. Because of our 
assumption on ms, for a large enough p c the field must 
start off very close to the minimum <j> = <f> c inside the 
body where the driving force V e g^ vanishes. Otherwise 
the solution (~ e msr /r) grows too fast inside the body 
and overshoots the exterior solution. For a sufficiently 
large p c the field stays close to <p c almost all the way to 6 
r = R. It follows from a second order Taylor expansion 
around 4>s that this is guaranteed to be the case as long 
as 



V c «,4>4>4>( ( t ) s,Pc){<i) o - (j>s) 
V e ff t <t,<t,((f)s,pc)msR 



(59) 



When all these conditions are satisfied, there exists of a 
critical solution in the limit p c — > oo which reads 

<\> = 4>c r<R, (60) 

{<t>c. - ^>oo)i? 



r> R, 



(61) 



which, apart from the numerical value of cfroo and (j> c , 
is completely model independent. This critical solution 
and its implications, for the case of power-law chameleon 
theories, was discussed in [9]. Another regime which 
can be described by exact solutions without having to 
solve model dependent equations is realised when 0oo 3> 
/3oo'7ipi < f ) Ar. In this regime the theory is effectively linear 
and the solution reads 



6mpi \R 2 
3mpir 



3 ) r<R, (62) 

r> R, (63) 



where /3oo = /3(0oo). This is the same type of solution as 
found in Newtonian gravity and the fifth-force-to-gravity 
ratio on a test mass outside the body is 



It 
F G 



20, 



2 

oo' 



while for the critical solution we find 

F 4> _ 2(j 2 ( \<t>oo - 0c| 

F G Po ° ^^mpi^iv 



(64) 



(65) 



Comparing the two cases we see that the critical solution 
corresponds to a screened fifth-force given that 



< 2/3 0O mpi$ 



A'. 



(66) 



which is exactly the screening condition we have found 
for chameleons, symmetrons and dilatons by solving the 
field equation explicitly. It is easy to show that the as- 
sumptions we started with do hold for these models. The 
critical solution, which formally only holds in the limit 
p c — > oo, will be a good approximation for the case of 
finite p c as long as the screening condition holds by a 
good margin. As current local gravity experiments give 
very tight constraints, if one wants to have cosmological 
signatures i.e. (3^, = 0(1), then this will be true in most 
cases. 

For the case where l^oo — </> c | ~ 2mpi/3 oa <fr]\r we would 
have to solve the model dependent equation to get accu- 
rate solutions. These solutions will interpolate between 
the two regimes found above, see e.g. [57] for a thorough 
derivation of chameleon equations in all possible regimes. 



C. Cosmological Scalar Field Dynamics 



5 As p c — > oo we have <f>s — > <f>c] the minimum for the matter 
density p c . The reason we explicitly write the limit here in- 
stead of taking <j>g = <f> c directly is to account for models where 
Mm^^^Ves,^ = °. but where lim Pc ^ 00 V o ff,<^(0s(Pc), Pc) = 
cxi as can be the case for generalised symmctron models as we 
shall see later on. Loosely speaking we can state this condition 
as: the mass at the minimum inside the body is increasing with 

Pc- 

For chameleons the solution only grows in a thin-shell close to 
the surface, but for large enough densities the field hardly moves 

at all. 



Here we consider the cosmological evolution of the 
scalar field (f> in modified gravity models with a mini- 
mum of V e s{<j>) at which the scalar field mass m satisfies 
m 2 ^> H 2 . The cosmology of the scalar field is tightly 
constrained by BBN physics due to the coupling of the 
scalar field to matter particles. The fact that the scalar 
field evolves along the minimum of V e e(cf>) implies that 
the masses of fundamental particles 



rriip = A(4>)m hat 



(67) 



in which mbare is the bare mass appearing in the matter 
Lagrangian, evolve too. In practice, tight constraints on 
the time variation of masses since the time of BBN 






p 



mpi' 



(68) 



where A<j> is the total variation of the field since BBN, 
impose that Am^/m^, must be less than ~ 10%. At a 
redshift of order z e sa 10 9 , electrons decouple and give a 
"kick" [10] to the scalar field which would lead to a large 
violation of the BBN bound. To avoid this, the field must 
be close to the minimum of V e ${4>) before z e and simply 
follow the time evolution of the minimum given by 



dV 



I Pm, 
TOpi' 



(69) 



Moreover, the total excursion of the scalar field following 
the minimum must be small enough. In practice, we will 
always assume that |</>/mpi| <C 1 along the minimum 
trajectory, implying that the BBN bound for the time 
dependent minimum is always satisfied. The models are 
then valid provided the electron "kick" does not perturb 
the minimum too much. We analyse this now. 

The background evolution of the scalar field is gov- 
erned by the homogeneous scalar field equation 



3^+^ = 0. 



(70) 



We assume that the contribution of the scalar field to the 
Hubble rate in the Friedmann equation is negligible until 
the acceleration of the Universe sets in 



where 



H 1 



1 ■ 



Prad + Pm + , 

3m P1 



= -^ + [A( ( f>)-l}p m + V (</>). 



(71) 



(72) 



The models that we consider here have a dynamical min- 
imum located at 4> m i n (t) such that 



dV cS : 



= 0. 



(73) 



Defining 5cf> = <fi — <j) m i n , we have for linear perturbations 
around the minimum 



6(f> + 3HSc/) + m 2 5<j) = F, 



(74) 



and the forcing term is then 

3p m oa~ 3 d 



F 



mpi dt 



A/3H 



(77) 



We must also take into account the "kicks" that the 
field receives every time a relativistic species decouples. 
These "kicks" correspond to the abrupt variation of the 
trace of the energy momentum tensor of a decoupling 
species at the transition between the relativistic and non- 
relativistic regimes. The abrupt change of T*} for the de- 
coupling species happens on a time scale much smaller 
than one Hubble time and can be modelled out using an 
"instantaneous kick" approximation [10] where the con- 
tribution to the scalar field equation is a (^-function. For 
kicks at the decoupling times tj , the source term becomes 



F = — 



3po d 
ropja 3 dt 



A/3H 



Aj3 ^ KjHjmpiStt -tj), 



(78) 
where Kj ~ <?«/<?* (toj) < 1 depends on the number of 
relativistic species g ir (rrij) at time tj and the number of 
degrees of freedom of the decoupling species gj . 

Let us now go through the different cosmological eras. 
During inflation, the Hubble rate is nearly constant and 
the field is nearly constant 7 . Indeed, the trace of the en- 
ergy momentum tensor is 



T« -12H 2 m 2 



pi, (79) 

in which p m = — p m = ZH 2 riip X is nearly constant in the 
slow roll approximation. As a result, the source term in 
the perturbed scalar field equation vanishes, and averag- 
ing over the oscillations with the fast period 1/m ^C 1/ H 
we have 

{5(j) 2 ) ex a" 3 , (80) 

implying that the field reaches the minimum of the effec- 
tive potential very rapidly during inflation. 

Assuming that reheating is instantaneous and that the 
field is not displaced during reheating, the field starts 
in the radiation era at the minimum of the effective po- 
tential during inflation. As the minimum has moved to 
larger values, the field rolls down towards the new mini- 
mum, overshooting and then stopping at a value 



^overshoot 



flation + i/6filmpi, 



(81) 



depending on the initial density fraction VL\ in the scalar 
field [10]. After this the field is in an undershoot situation 
where the field is essentially moved according to the kicks 



where 



3Hcj> = -Af3 V KjHjmpiS (t -tj), 



(82) 



F 



l__d 
a?~dl 



dt 



Using the minimum equation, we find that 

0min = ?PA , 

m z mpi 



(75) 



(76) 



7 Note the paramctcrisation m(a) = moa~ r to be introduced be- 
low only applies when the scalar field is sourced by the pressure- 
less matter, and does not apply to the inflationary era, in which 
<f> remains nearly constant simply because the density of the in- 
flaton docs so. 



Each kick brings the field to smaller values, with a vari- 
ation 



A& 



-ftjAjKjmpi, 



(83) 



in the radiation era [10]. Although the details depend on 
the kicks and the initial energy density of the field, we 
can assume that after all the kicks before BBN, the field 
is close to the minimum of V e g ((/))■ We will assume that 
this is the case by z; n i « 10 10 where the matter density 
is equivalent to the one in dense bodies on Earth today. 
If this were not the case then the field would move by 



A0e 



-PeAeKeTTlpi, 



(84) 



when the electron decouples during BBN, and the masses 
of particles would vary too much during BBN. Note that 
for the rest of this subsection a subscript e will be used to 
denote the value of a quantity at the electron decoupling. 
Hence viable models must be such that the scalar field 
remains in the neighbourhood of the minimum since well 
before BBN. In this case, the deviation of the field from 
the minimum can be easily obtained from 



3H5(f> 



m 5<fi 



3pmO a 



TOPI 



(If 



A/3H 



-A e /3 e K e H e mpi5 (t - t e ) , (85) 



where we only take into account the electron kick. Dcfin 
ing 8(j) — a -3 / 2 ?/;, we find that 

3p m0 a- 3/2 d \A0H' 



v> 



9w h2 
4 



i> = - 



mpi di 
A e j3 e K e H e a 3 J 2 m P i5 (t - t e ) 



(86) 



As m 2 ^> H 2 , the solution is obtained using the WKB 
approximation and reads 



APH 



5(f) 9O m0 ifg d 

mp\ a 3 m 2 di 



-6 (t - t e ) A e P e K ( 



H, 



3/2 



y/m e m a 3 / 2 



(87) 



m(t')dt', 



in which the second term is only present when t > t e , 
being the Heaviside function. We will always assume that 
j3 and m vary over cosmological times, hence we have 






A/3H 



»? z 



m 



ApH 2 



(88) 



in which g(t) is a slowly- varying function of time whose 
value is of order unity. Averaging over the rapid oscilla- 
tions, we have 



SIQ^A*? H* mt H* 



'pi 



772n TO TW 



AlBlnl al H 2 m e 



ei ' e e ^e 



(89) 



The first terms is of order ft^H^/m^ <C 1 now, implying 
that it has a negligible influence on the particle masses. 



This guarantees that the minimum is indeed a solution 
of the equations of motion. The second term corresponds 
to the response of the scalar field to a kick. It is initially 
very small as suppressed by H 2 /m 2 <ti 1, implying a tiny 
variation of the fermion masses during BBN. Its influence 
increases with time as 1/ma 3 and we must impose that 
this never compensates the fact that H 2 /m 2 is extremely 
small. 

Consider an interesting example with m{a) = moa~ r 
which will reappear later. In such a case the second term 
in the above equation can be rewritten as 



A 2 J 2 e n 2 e alH 2 m e 



Uq il 



TTIq Sl m o 



r0 or-V- 3 , 



(90) 



where we have assumed A 2 fi 2 n 2 ~ 0(1) and il r o <C fi m o 
is the fractional energy density for radiation (photons and 
masslcss neutrinos) at present. From this formula we can 
easily see that 

1. when r < 3 the minimum of V e s given by the min- 
imum equation is an attractor, because the magni- 
tude of the oscillation decreases in time; 

2. assuming that Hq ~ 10~ 3 mo (see below) and 
U m o ~ 10 3 ri r o, then today we have (S<j) 2 )/m Pl ~ 
10 _9 ae _1 which is of order one if r — 0. Clearly, for 
r < 2 the amplitude of oscillation can be too big 

{\/(S(j) 2 ) > </>min) at early times; 

3. if r > 3 which is the case for f(R) gravity models in 
which f(R) ~ R + Rq - Rx{R tr /R) n , s/T5^)/mp\ 
increases with time but never becomes significantly 
large. For example, if r = 3 then y/(8<f> 2 )/mpi ~ 
10 _15 today, which means that, although the mini- 
mum of V e ft(4>) is not strictly speaking an attractor, 
it is extremely stable to kicks and governs the back- 
ground dynamics of the model. 



D. The Equation of State 

We have described how the cosmological constraint 
from BBN imposes that the scalar field must be at the 
minimum of the effective potential since BBN. As such 
the minimum of the effective potential acts as a slowly 
varying cosmological constant. We have also seen that 
when m 2 ^S> H 2 , a large class of models are such that 
the minimum is stable. In this case, the dynamics are 
completely determined by the minimum equation 



dV" 



-f3A 



Pm 

mpi' 



(91) 



In fact, the knowledge of the time evolution of the mass 
m and the coupling )3 is enough to determine the time 
evolution of the field. Indeed, the mass at the minimum 
of V e s, 



m 



d 2 y eff w 



(92) 
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and the minimum relation leads to 

a4> z m P1 d(p mpi 

where the couplings to matter /3 can be field dependent. 
Using the minimum equation, we deduce that the field 
evolves according to 



d(/) _'SH p m 

17 ~ — 2" • 

dt m z mpi 



(94) 



This is the time evolution of the scalar field at the back- 
ground level since the instant when the field starts being 
at the minimum of the effective potential. In particular, 

we have 



1 /# 

2 V di 



2 = 2 la m p 2 A 2 (-X P „ 

2 \m 1 



(95) 



which is tiny compared to p m . 

Because of the interaction between the scalar field and 
matter, the energy momentum tensor of the scalar field 
is not conserved. Only the total energy momentum 

ptot = -3H(ptot +Ptot) (96) 

is conserved, where the total energy density is 

Ptot = Pm + P<j> (97) 

with 



p<p = y + v cS {(f>) 7 

ft 
ptot = P4> = y -v{(j>), 



(98) 
(99) 



and where we have neglected the radiation component 
in the matter era. It is crucial to notice that the energy 
density of the scalar field involves the effective potential 
V^ff while the pressure only involves V. This is a crucial 
feature of scalar-tensor theories. 

We can define the effective equation of state of the dark 
energy fluid as 



Wd, = — ■ 
P<t> 



(100) 



Using the Friedmann equation we find the Raychaudhuri 
equation involving the effective equation of state w<f, as 



1 



a 6rnpj 



6to p1 



[Pm + (l + 3u>0)p ] 
(1 + 3u>tot)/9tot 



where we have defined the total equation of state 

Ptot 



Wtot 



Ptot 



(101) 



(102) 



The universe is accelerating provided a > which leads 
to 



1 

Wtot < —Z 

as expected, which is equivalent to 

«V< -- 1+ — 
3 V P4> 



(103) 



(104) 



The situation of the modified gravity models can be easily 
analysed as 



■w<t> + 1 



{A-l) P r, 



f + v{4>) + (A- IK 



which can approximated as 



ft , , s a 

w + 1 « T77-T + (A - 1) 



v{4>) 



a 



(105) 



(106) 



The first term corresponds to the usual quintessence con- 
tribution and the second term can be approximated as: 
1±^ 2__£*_<W = 3p2Q h^o^ This j Ues 

mpi ilj, mpi V.^4, U4, ^ m* SZ^ ' 

that 

W(j> + 1 „ (A- 1)^ « 3tt m /? 2 f-V 7^- (107) 

In the recent past of the Universe where Q m and O^ are of 
the same order of magnitude, this implies that the back- 
ground scalar field acts as a cosmological constant due 
to the large H 2 /m 2 suppression. In the past, the back- 
ground cosmology deviates from a ACDM model only if 
fi^ becomes so small that it compensates m 2 /H 2 . We 
will not consider this situation in the following. 



III. MODIFIED GRAVITY TOMOGRAPHY 

A. Reconstruction of the Dynamics 

We have seen that when m 2 ^> H 2 a large class of mod- 
els are such that the minimum of the effective potential 
is stable or quasi-stable, and in these cases the dynamics 
are completely determined by the minimum equation 



dU 



-0A 



Pm 

mpi' 



(108) 



In fact, the knowledge of the time evolution of the mass 
m and the coupling j3 is enough to determine the bare po- 
tential V(4>) and the coupling function A(<fi) completely. 
To see this, integrating Eq. (94) once, we find 



<Ka) 



TTipi Jo, i am 2 (a) 



p m (a)da + 4> c , 



(109) 



where (f> c is the initial value of the scalar field at aim < 
a BBN an d we have taken A(<f>) wl as the temporal varia- 
tion of fcrmion masses must be very weak. If the coupling 
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/3 is expressed in terms of the field 
factor a, this is also equivalent to 



and not the scale 



m 



3 
mpi 



1 



am 2 (a) 



p m (a)da. 



(110) 



Similarly the minimum equation implies that the poten- 
tial can be reconstructed as a function of time 



V = V n - 



P 2 (a) 

m Pl Jaini am2 ( a ) 



p 2 m (a)da, 



(111) 



where Vq is the initial value of the potential at a = a ln \. 
This defines the bare scalar field potential V(0) para- 
metrically when /3(a) and 771(a) are given. Hence we have 
found that the full non-linear dynamics of the theory can 
be recovered from the knowledge of the time evolutions 
of the mass and the coupling to matter since before BBN. 



B. Tomography 

The previous reconstruction mapping gives a one-to- 
one correspondence between the scale factor a and the 
value of the field 0(a) in the cosmic background. As the 
scale factor is in a one-to-one correspondence with the 
matter energy density p m (a), we have obtained a map- 
ping p m —J- 4>{p m ) defined using the time evolution of 
m(a) and /3(a) only. Given these evolutions, one can re- 
construct the dynamics of the scalar field for densities 
ranging from cosmological to solar system values using 
Eq. (109) and Eq. (111). By the same token, the inter- 
action potential can be reconstructed for all values of 
(and pm) of interest, from the solar system and Earth 
to the cosmological background now: a tomography of 
modified gravity. 

In particular, we can now state the screening condition 
of modified gravity models as 



/?(«) 
am 2 (a) 



p m {a) Aa < /3 utro P1 $jv, 



(112) 



with constant matter densities Pm,out = Pm{o, = ai n ,out) 
inside and outside the body respectively, and where we 
have defined /3 out = /3(a — a out ). It is remarkable that the 
gravitational properties of the screened models are cap- 
tured by the cosmological mass and coupling functions 
only. 



C. Dilatons 

Let us consider a first example: the dilaton models in 
which the coupling function /3(0) vanishes for a certain 
value 0* of the scalar field 0. On the other hand, we 
assume that the potential is positive definite and is of 
runaway type. It is enough to study the dynamics in the 
vicinity of the field 0» , where 



/3(0) sa A 2 mp\((f)- 0* 



(113) 



from which we deduce that 

In 



= 9A 2 ml i n m0 H$ / 



da 



and therefore 



|/3(0)| = \0(<j> c ) I exp 



9A 2 m,p l Q [ 



lOHn 



a 4 m 2 (a) ' 



da 



(114) 



a^m 



'(a) 



(115) 
In particular, we find the relation between the coupling 
at the initial time and other cosmological times. 

The initial coupling (taken at a m i < obbn) is the same 
as in dense matter on Earth, as long as the field minimises 
its effective potential in a dense environment, and it is 
related to the cosmological value of /3 today, /3(0o), by 



|/3(0o)| = |/3(0 c )|exp 



9A 2 ml i n r 



X)H Q 



da 



a 4 m 2 (a) 



(116) 

It is possible to have a very small coupling in dense mat- 
ter |/3(0 C )| <C 1 for any value of the coupling on cosmo- 
logical scales |/3(0o)| provided that A 2 > and that the 
time variation of m(a) is slow and does not compensate 
the 1/a 4 divergence in the integrand. In this situation, 
the coupling function /3 converges exponentially fast to- 
wards zero: this is the Damour-Polyakov mechanism [8] . 
The fact that A 2 > guarantees that the minimum of 
the coupling function is stable and becomes the minimum 
of the effective potential which attracts the scalar field in 
the long time regime. If A 2 < 0, the effect of the coupling 
is destabilising and implies that diverges exponentially 
fast away from 0*. 

Alternatively, a smooth variation of the coupling func- 
tion to matter in the cosmological background and there- 
fore interesting consequences for the large-scale structure 
can be achieved when the evolution of the mass of the 
scalar field compensates the 1/a 4 factor in the radiation 
era and evolves in the matter era. This is obtained for 
models with 



m 2 (a) = 3A 2 H 2 (a)ml 



pi- 



(117) 



Indeed, H(a) 



in the radiation era, which implies 



that the time variation of /3 between BBN and matter- 
radiation equality is 



/3(0)=/3(0 c )exp 



O— (a- Oi„i) 



and in the matter dominated era 



>(0c q ) — ) =/3(0e q )- PmKq) 



Pra{a) 



(118) 



(119) 



where a subscript e q denotes the value of a quantity at 
the matter-radiation equality This is the behaviour of the 
dilaton models we have already analysed gravitationally 
in 8 IIB2. 
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D. Symmetron 



E. Generalised Symmetrons 



In the symmetron models the coupling to matter van- 
ishes identically in dense regions or at redshifts z > z*, 
while a larger coupling is obtained after a transition at 
a redshift z* and in the low matter-density regions. This 
can be obtained by choosing 

(3(a) =P^ I -{^) 3 , (120) 

for z < z+ and /3 = 0, z > z+. Similarly we choose 



m(a) = m±\l 1 — ( — J . 



(121) 



Using the reconstruction mapping, it is straightforward 
to find that 



m-W-(t 



(122) 



for z < z± and cf> — before. The potential for z < z+ 
as a function of a can then be reconstructed, using the 
technique introduced above, as 



V(a) = V + 



2m2mpj 



(t)' 



where 



PmO 



(123) 



(124) 



is the matter density at the transition between </>(a) = 
and 4>(a) > 0. The potential as a function of <f> is then 



V{<j>) = v + V - ^-</> 2 , 



where 



and 



together with 






l* = V2/i, A = — , 



m = ^ 



(125) 



(126) 



(127) 



(128) 



This completes the reconstruction of the particular sym- 
metron model presented in [17] from m(a) and /3(a). 



With the parametrisation developed in this paper it is 
easy to create new models (in a more intuitive way than 
starting with the Lagrangian) by changing the mass and 
coupling functions. Here we give a simple example by 
generalising the symmetron models. 

We start by generalising the coupling function 
Eq. (120) 



/3(a) = & 



(1 



a*\ 

a 



1/8 



(129) 



for z < z± and /3 = for z > z±. Similarly we choose 

i/p 



m(a) = m* 



where the field evolves as 



1 



(t) 



4>(a) = </)* 

where we have defined 

2(p-q+pq) 

m = — -. 

p-2q + pq 



-I — 

a 



2p-2q + pq 
p-2q + pq ' 



and where 



(m - n)/3*P* 



mjmpi 



(130) 



(131) 



(132) 



(133) 



Eventually we find 

(m - n)Plpl 



V{<t>) = V + 



and 



'*'"pi 



P(4>) = & 



1 

m \4>* 



(134) 
(135) 



The indices to and n should be taken to be even integers 
to keep the potential symmetric around <fr = 0. The stan- 
dard symmetron corresponds to the choice to/2 = n = 2. 
We can now show explicitly that this generalised sym- 
metron model has the screening property as we did for 
the original symmetron model in § IIB2. Let us consider 
a spherically dense body of density p c and radius R cm- 
bedded in a homogeneous background. The field profile 
inside the body is 



4> = 4>s 



sinh TOgr 
msr 



r < R 



(136) 



where 



V d(j> J s mpi 



n — lp c 



n p+ \4> 



n-2 



^ ( 137 ) 
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is the scalar field mass at r = 0, <j>s the corresponding 
field value and p* is as in the symmetron model the crit- 
ical matter density when the transition of the minimum 
of V e e{4>) from (f> = to <fr — ±</>* takes place in the 
cosmological background. 

The field outside the body, on scales shorter than the 
large range m^ 1 , is 



Local Constraints 



D 
r 



r>R 



Matching at r = R gives us the solution 

4> s cosh(ms-R) = 4>* 

/tanh(TOsi?) 



D 



KR 



(138) 

(139) 
(140) 



\ msR 

The first condition, which determines (f>s, can be written 

n/2— 1" 



9S , 
-— cosh 



1 



(141) 



where a = ? "_^ y- (m+R) 2 . We can change it into a simple 
equation for msR 



im s R) 2 cosh™" 2 (m s R) = a 



(142) 



From these equations we see that when a » 1 we get 
4>s ~ 0, msR ^> 1 and therefore D ss —<fi*R. Note that 
if n > 2 the mass vanishes at <\> — 0, however, this is not 
a problem for the screening mechanism. Even though a 
large a pushes the field down towards — 0, tos is still 
an increasing function of a according to Eq. (142). 

The fifth-force on a test mass outside the body is found 
to be screened as long as 



|0c-</>oo| < 2m P i/3*$ 



N 



(143) 



where <j) c = 4>s ~ and 0oo = 0*. This condition is equiv- 
alent toa> 1 and shows that the screening property is 
present in this model. 

Comparing the case n = 2 with n > 2 we find that even 
though (f>s/<P* is larger in the latter case, the coupling 
f3{4>s) is smaller as long as we have screening. This means 
that the force between two test-masses in a dense envi- 
ronment is more screened for larger n. Local constraints 
for the generalised symmetrons are therefore satisfied for 
(at least) the same range as the standard symmetron: 
m /H > 10 3 . 



IV. RECONSTRUCTING f(R) MODELS 

A. Gravity Tests and Chameleons 

Consider now the important case of a non-vanishing 
coupling function /3(a). Defining /3(a) = flog (a) and m = 
toq /(a), we find that 



TOpi 



9/3 n 



mO" 



Hi 



(1(7 



9 ( a ) 
a 4 / 2 («) ; 



(144) 



10 s - 



1 o"- 



110 




FIG. 1. The constraints on mr, /Ho as a function of r for /3o = 
l/\/6 and s = 0. Valid models must be above the red (solar 
system), mauve (cavity) green (m > H), light blue (mL > 1), 
light red (fi) and brown (galaxy) lines. The blue line gives the 
detectability of effects on the CMB by the Planck satellite. 
The strongest constraints are the cavity and galactic bounds 
for small and large r respectively. Models with r > 3 satisfy 
the constraints and can lead to a modified gravity regime on 
large scales. 



which allows one to test the screening properties of these 
models. 

Let us first consider the solar system tests. Evaluating 
Eq. (144) in the galactic background, we find that 8 



TOpi 



= 9/3ofi ro o^f r\\a- 9{a) 



'0 Ja iai 



a 4 / 2 («) : 



(145) 



where a ga i w 10~ 2 is the scale factor when the matter 
density in the cosmological background equals the galac- 
tic density p ga i ss 10 6 p c . Defining 



Ogal 



AR 



R 6to p1 £ c $ g 



(146) 



where R is the radius of a spherical body, the modifica- 
tion of gravity in the solar system has a strength 



2/3gal/3 c 



3AR r . 



R. 



(147) 



In this expression /3 ga i is the value of the coupling func- 
tion /3(4>) in the galactic background, $© is the value of 



Again, here for simplicity we have assumed that the scalar field 
minimises V e g((j)) in the galactic background. While this is true 
for certain parameter space, in general it should be tested against 
numerical simulations. 
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the Solar Newtonian potential ($ Q ~ 10 -6 ) and /3 C is the 
coupling inside a dense body. The magnitude should be 
less than 10 -5 to comply with the Cassini bound in the 
solar system [50] . This condition is independent of j3 c and 
reads 



m 



gal 



The integral 



da 



7(a) 



a 4 / 2 (a) 



<10" 



9Cl m oH 



;$c. 



1 = 



da- 



9W 



(148) 



(149) 



a 4 / 2 (a)' 

is potentially divergent for small values of a; n ; ~ 10 -10 . 
Hence we must impose that f{a) 2 /g(a) compensates the 
1/a 4 divergence in the integrand. As mentioned above, we 
have assumed that galaxies are screened to minimise the 
disruption of their dynamics, although the necessity of 
this condition should be ascertained using iV-body sim- 
ulations [31]. Enforcing the screening condition imposes 



l^gal - <ft)| < 6/3o"1Pl$gal, 



(150) 



in which the galactic Newtonian potential is $ ga i ~ 10 6 
and 



^0 - ^gal _ QR „ Ho 

mpi mf 



da- 



.9 (a) 



J Ogal 



« 4 / 2 («)' 



(151) 



A slightly stronger bound is obtained from the Lunar 
Ranging experiment [49] with the 10~ 5 on the right-hand 
side of Eq. (148) replaced by 10~ 7 . 

Strong constraints can also be obtained from labora- 
tory experiments. Using the fact that the initial matter 
density at z; n i ~ 10 10 is roughly the same as that in a typ- 
ical test mass in the laboratory, gravity is not modified 
provided test bodies are screened, i.e., 

|</>lab - <t>c\ < 2/3 c TO P i$ lab , (152) 

where $i a b ~ 10~ 27 for typical test bodies in cavity ex- 
periments of size -L, and </>i a b = (f>{ai a b) is determined by 
"i(aiab) ~ 1/L (see the appendix for more details). 



B. f(R) Gravity Reconstruction 

Viable f{R) models are nothing but chameleons [32] 
with a constant value of the coupling function fi{4>) = 
1/V6- We have already described the background dy- 
namics of these models. Here we shall derive the mapping 
between the evolution of the scalar field mass m(a) and 
the function f(R) for curvature values ranging from the 
ones in dense bodies to cosmological ones. These models 
are equivalent to chameleon models where the potential 
is given by 9 



V{(j>) = m P1 - 



2/1 



(153) 



9 In the discussion of f(R) gravity we shall use R to denote the 
Ricci scalar. 



in which fn = df /dR. The mapping between R and <f> is 
given by 



f R = cxp -2/3 



-^ 



mpij 

Given the mass function m(a), we have 

da 



and 



(a) = 9(3n, m0 H m Pl 



V = V Q 



„ ini a A m?{a) 
P 2 P 2 m {a) 



da. 



/ Q . n . am 2 (a) m P1 
We can reconstruct R(a) using the fact that 



R{<f>) 

and f(R) using 



2/3 m 

e m pi 



i_ 1 



/3mp\ d(f> 



2rf- 



e~ Ap ^V{4>) 



(154) 



(155) 



(156) 



(157) 



/(i?) = i?(0)e- zp ^T- — 



-4/3- 



~V(4>), (158) 



'pi 



which is equivalent to 



f(R) = 



-4/3- 



■v(4>)- 



-4/3 



J^dV 



'PI 



/?mpi 



(159) 



once we have obtained V{<f>) from the above implicit pa- 
rametcrisation. 

When (3<fi/mp\ <C 1 as required from the BBN con- 
straints, the above equations can be simplified and read 



f(R)=R-2^ 



m 



Pi 



where 



1 (JU 4 

R{<t>) = — R — xr + — W 

pmpi d<p rnpj 



(160) 



(161) 



This is the parametric reconstruction mapping of f(R) 
models. 



C. Large Curvature f(R) Models 

We can apply these results to the case with m = rriQa~ r 
leading to models where 



9Q m o/3H 2r _ i 



mpi (2r-3)m 2 im 
which reduces to 



a 

aini 



2r-3 



9n m0 (3H 2 



mpi 



(2r - 3)mg 



o a 2r-3 



, (162) 



(163) 
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at late times. Similarly we have 
V(a) = V - 



W 2 P 2 



>nO 



,2r-6 „2r-6 



2(r - 3)mpjm5 
Now for late enough times we have 

V = Vo - C 



";, 



mpi 



2(r- 3) 
2r-3 



(164) 



(165) 



for a constant C . Notice that for 3/2 < r < 3, these mod- 
els are chameleons with an inverse power law potential 
V{<j>) ~ 0"™ with 



„ 3-r 

n = 2 

3-2r 

We can equivalcntly find that 

m 2C r - 3 " 



13ml, 2r ~ 3 



Finally we find that 



mpi 



Vb 



(166) 



(167) 



'pi 



/(;?) = R _ 

m P1 



K, + C 



'i?-4 



U 



I?* 



(168) 



where i?* = 2(r - 3)C/ [(2r - 3)/3m P1 ] and 



n= g( r -3)- 



(169) 



Large curvature models are defined for r > 3 here. This 
completes, in this particular example, the reconstruction 
of the f(R) models from the knowledge of the function 
m(a). 

The gravitational constraints for these models have 
been fully analysed in [5]. We have summarised these 
constraints in Figure 1 where we see that the strongest 
constraints on the range of the scalar interaction arise for 
r < 3, i.e., for inverse power law chameleon models. For 
r > 3, i.e., for large curvature f(R) models, the screening 
of the Milky Way is a loose constraint which needs to be 
further analysed with TV-body simulations. 



allowing one to reconstruct the field history entirely: 



In. - fm = 18/3 fl m oH Q 



1 



a a 4 m 2 (a) 



da, (172) 



which depends on the mass evolution uniquely. This can 
be rewritten using the i?-function. In fact, using 



d# 3,„ , rri 

— = --(l + w)Hdt, 



(173) 



in an era dominated by a fluid of equation of state w, we 
find that 



B 



Irr I R 

J R 2,{l + w)H' 



23- 



With fn = e ™pi we have 



AR f R d0 

JRR-TT = ~ l P "77 

at mpi di 



and therefore 



B = 



AP dej) 



3(1 + w)m P i Hdt' 

and using the minimum equation we get 

B = — — Vl m — -. 

1 + w m^ 



(174) 



(175) 



(176) 



(177) 



Because f3 = l/v6, in the matter dominated era this 
gives 



B — Q, m - 



H 2 _ 

2 



in 



(178) 



which is completely determined by m(a). Hence we find 
that 



Jr - Jro = 3 / — da. 



(179) 



The knowledge of B(a) and Jrq determines the back- 
ground evolution in the f(R) gravity models in a com- 
pletely equivalent way to the m(a) parameterisation. 



V. GROWTH OF LARGE-SCALE STRUCTURE 



D. Comparison with the B-parameterisation 



The f(R) theories arc generally parameterised using 



[29] 



B=l^H^. 
Jr dH' 



and Jr — 1 now. As 0/mpi ^C 1 we have that 



f R - 1 = -2/3- 



mpi 



(170) 



(171) 



We have shown that the non-linear structure of the 
screened models can be reconstructed from the knowl- 
edge of the mass and coupling functions. These functions 
are time dependent only. In particular, we have seen that 
this allows one to fully analysed the gravitational tests 
and the cosmological background evolution. Moreover we 
have shown that the cosmological dynamics typically is 
indistinguishable from a ACDM model at the background 
level. Here we will find that this is not the case at the 
perturbative level and that the mass and coupling func- 
tion allow a full description of the linear and non-linear 
regimes. 
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FIG. 2. The relative difference of the matter power spectrum P(k) in the chameleon model from that in the ACDM model with 
exactly the same background expansion history, initial conditions and physical parameters. Upper Left panel: the dependence 
of the result on the modified gravity parameter j3o. Upper Right panel: the dependence of the result on the parameter r. Lower 
Left panel: the dependence of the result on the parameter s. Lower Right panel: the dependence of the result on the parameter 
£ = H /m . 



Linear Structure Growth 



The linear perturbation equations for a scalar field cou- 
pled to matter particles are listed in [56] in the covariant 
and gauge invariant formalism. Denoting by A m the den- 
sity contrast of the pressure-less matter, v m its velocity 
and S(f> the perturbation 10 in the scalar field, their evo- 



10 Note that this is different from above, where we used S<f> to denote 
the oscillation of the background <j> around m i n (£) . 



lution equations are as follows: 



1 Pn 



A" +_ A' - -i2Ii_ a 2/\ 
a Z ?rip| 

+k/3{a)mpl {k5(jy - <p'v m ) = 0, (180) 
v' m + —v m + /3(a)m p { (<f>'v m - kS(f>) = 0, (181) 



a m 



'(a)] 6<f> 



-/3(a)^-a 2 A m + k(t>'Z 
mpi 



0, (182) 



where a prime denotes the derivative with respect to the 
conformal time, kZ = \t' in the Newtonian gauge is a 
variable of the curvature perturbation which is irrelevant 
for our discussion since it is multiplied by c6'/mpi <ri = 
a' /a, and we have neglected contribution from radiation 
as we are focusing on late times. 

Neglecting the terms proportional to (f> in the above 
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FIG. 3. The relative difference of the matter power spectrum P(k) in generalised symmetron models from that in the ACDM 
model with exactly the same background expansion history, initial conditions and physical parameters. Upper Left panel: 
the dependence of the result on the parameter a* - the scale factor value at which the symmetry breaking of the effective 
potential happens. Upper Right panel: the dependence of the result on the modified gravity parameter /3*. Lower Left panel: the 
dependence of the result on the parameter q. Lower Right panel: the dependence of the result on the parameter £ = Ho/m+. 
As an example we have chosen p = 2. 



equations we get the following equation [10] 



A" 



-A. 



1 Pm 

2 m|j 



o/A r 



2f3 2 (a) 



IM 



= 0, 

(183) 

where we have used the fact that, given that in Eq. (182) 
the term k 2 + a 2 m 2 3> H 2 , 6~<f> follows the solution 



P{a) 



Pn 



k 2 + a 2 m 2 (a) mpi 



a 2 A r 



(184) 



and rapidly oscillates around it (see more details below) . 
On very large scales, k <C am(a), we can see that 
Eq. (183) reduces to 



A" 



1 Pn 



-a 2 A r 



0. 



(185) 



'pi 



which governs the growth of matter density perturbation 
in the ACDM model. The effect of modified gravity is in- 



corporated in the second term in the brackets of Eq. (183) 
and becomes significant when am(a)/k < 1, namely for a 
light scalar field mass 771(a) or on small length scales. For 
all models shown here the CMB spectrum is the same 
as the ACDM prediction, because the scales relevant for 
the CMB is very large and therefore not affected by the 
modified gravity. 

In order to illustrate these considerations, we have 
computed the linear matter power spectra P(k) for a 
number of generalised chameleon (Fig. 2) and symmetron 
(Fig. 3) models. 

For the generalised chameleon models, we have used 



m a r , (3 = f3 a~ 



(186) 



The impact of gravity tests for (3 = 1/V&, s = have 
been given in Fig. 1. There we can see that values of r > 3 
are favoured by the local gravity tests. We have varied 
the four parameters in the paramcterisation of /3(a) and 
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m(a): /3o,r,s and too- Because mo is not dimensionlcss, 
we have defined a new variable £ = Hq/ttiq instead. We 
find the following results, all as expected: 

1. increasing the coupling /3 strengthens the modifi- 
cation of gravity, which causes more matter cluster- 
ing, resulting in a higher matter power spectrum; 

2. r characterises how fast the scalar field mass de- 
creases in time: the higher r the faster it decays. 
Given that too is fixed, a higher value of r means 
that the Compton wavelength (essentially the range 
of the modification to gravity) decreases faster in 
the past, and therefore the modification of gravity 
starts to take effect later - this would mean less 
matter clustering; 

3. s specifies how fast the coupling function changes 
in time: s = implies /3(a) remains constant, while 
s > (s < 0) means fi(a) decreases (increases) 
in time. If (3q is fixed, the larger s is, the larger 
/3(a) becomes at high redshifts - this would mean a 
stronger modification to gravity and stronger mat- 
ter clustering; 

4. £ specifies how heavy the scalar field is, or equiv- 
alcntly the range of the modification of gravity: 
smaller £ means shorter Compton length of the 
scalar field, and therefore weaker matter clustering. 

The potential of the generalised symmetron models 
has been given in Eq. (132,134), but one should be care- 
ful that the parameters p, q (or equivalent n, to) cannot 
take arbitrary values. For example, <f> n might not be well- 
defined if <j> < 0. Here let us consider the special case with 
p = 2 (n = 2, m = 2 + q), in which the potential becomes 



V(<f>) = v + ^$4- 



1 




2+q 


1 

" 2 


.<t>*. 


2" 
( 


2 + o 



'PI 

and this avoids the situation in which the scalar field 
becomes massless at <j> = 0. Furthermore, choosing q = 
2, 4, 6, • • • not only ensures that 4> 2+q is well-defined for 
any value of <f>, but also makes the potential symmetric 
about </> = 0, as in the original symmetron model. Finally, 
with p — 2 another property of the original symmetron 
model, that /3(<f)) oc </>, is preserved as well. 
Again, the results in Fig. 3 are as expected: 

1. increasing a* implies that the modification of grav- 
ity starts to take effect at a later time, and this will 
weaken the matter clustering; 

2. increasing /3* increases the coupling strength over- 
all, and leads to stronger matter clustering; 

3. increasing q increases /3(a) for a > a* and causes 
stronger structure growth; 

4. decreasing £, as in the chameleon case, decreases 
the range of the modification of gravity, and there- 
fore leads to less matter clustering. 



k = 1.0 h/Mpc 




0.1 
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FIG. 4. An illustration of the time evolution of the scalar field 
perturbation 5(f). The black solid curve is the numerical solu- 
tion while the green dashed curve is the analytical approxima- 
tion given in Eq. (184). The results here are for k — 1/iMpc" 
but the qualitative feature remains for other values of k. The 
modified gravity parameters are shown beside the curves. 



Before we finish this subsection, let us come back to the 
evolution of the scalar field perturbation 8(f). As explained 
above, an analytic approximation to this can be obtained 
in Eq. (184). However, as for the background evolution, 
where <j> oscillates quickly around </> m in(£), we may expect 
that the true value of 8(p oscillates around the analytic 
solution as well. This is confirmed in Fig. 4. 

In the model shown in Fig. 4 we have chosen r = 3.0. 
Obviously, the larger r is, the larger the scalar field mass 
m(a) becomes at early times. A rapid decrease of m(a) 
would mean that the effective potential for 8(j) changes 
its steepness very quickly. Suppose the oscillation of 8<p 
has some initial kinetic energy, then as the effective po- 
tential becomes less steep the amplitude of the oscilla- 
tions increases since the kinetic energy does not disap- 
pear quickly. Consequently, if we increase r further we 
get even stronger oscillations and if, in contrast, we de- 
crease r then the oscillations become weaker. We have 
checked explicitly that for r = 1.0 there is essentially no 
oscillation. 

At late times Ho/mo = £ ~ 10~ 3 , which implies that 
the period of the oscillation is roughly 10~ 3 the Hubble 
time, and is much longer than the typical time scales for 
human observations. As a result, one cannot average 8<f> 
over several periods to get (8(f)) . Indeed, as the amplitude 
of oscillation in Fig. 4 is bigger than the analytic solution 
of 8(j) in Eq. (184), the value of 8(j> one observes at a given 
time is rather random and could be far from the one given 
in Eq. (184). This is the case for the f(R) gravity model 
in [31], where r = 4.5. 

Whilst this seems to be a problem, this is not really the 
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FIG. 5. The time evolution of ^(k, a) for a chosen value of 
k = 0.1/iMpc -1 as an illustration. The black solid is the full 
numerical solution, the red dashed curve is obtained using 
the numerical value of $jv using the analytical solution of 8(f) 
given in Eq. (184), while the blue solid curve is Eq. (196). The 
modified gravity parameters are shown beside the curves. 



case. Indeed in the solar system the matter density is so 
high that the oscillation is faster than it is in the cosmo- 
logical background, and we actually observe the averaged 
value (6<fr) . On linear scales, as 8<j> oscillates, overshoot- 
ing and undershooting the value given in Eq. (184), we 
have checked by replacing the numerical solution of 5<j) 
by the analytical formula given in Eq. (184) that we ob- 
tain identical power spectra P(k) in the two approaches. 
Hence the mean value solution Eq. (184) gives a very 
good description of the statistical properties of linear per- 
turbations. 



B. The Jordan Frame Picture 

In this section we compare our results with a simple 
and effective way of paramctcrising linear perturbations 
which has been used in the literature in the past few 
years [38-45] (other interesting and more general ap- 
proaches for the linear regime include the parameterised 
Post-Friedmann framework of [36, 37] and the fully co- 
variant parameterisation of [46-48]). Such a way of pa- 
rameterising any modification of gravity utilises two arbi- 
trary functions fi(k, a) and j(k, a) through the (modified) 
Poisson equation 



fc 2 \f/ = ATni(k,a)GNa 2 5pr, 



and the slip relation 



$ = 7(fe, a)*. 



(188) 



(189) 



Here Gn is the bare Newton constant, and ^ and $ are 
the two gravitational potentials in the Newtonian gauge: 



ds 2 



'(1 + 2*)d?7 2 + a 2 (l - 2$)dx 2 



(190) 



in which (77, x) are the conformal time and comoving co- 
ordinates. 

So far we have focused on the Einstein frame. In the 
Jordan frame as described by the line element above, 
the perturbative dynamics can be described using two 
Newtonian potentials where we have the relation 



d3* = A z {4>)As z 



(191) 



and ds 2 is the line element in the Einstein frame ex- 
pressed in the Newtonian gauge. Expanding in perturba- 
tion around a background value with A[<^>(£)] ~ 1, we can 
relate these two potentials to the Einstein frame Newton 
potential 



f = $ 



JV 



P- 



mpi 



$ = $ 



A' 



/3- 



TOPI 



(192) 



Hence we see that in the Jordan frame the two Newtonian 
potentials are not equal, a fact which can be interpreted 
as resulting from the existence of a non-anisotropic stress 
contribution coming from the scalar field. It is useful to 
define 



e(k,a) 



2f3 2 



1 



(193) 



k 2 



Using the definitions in Eq. (192), the analytical approx- 
imation for S(f> in Eq. (184) and the Poisson equation 



- k H N = 1 --^a 2 A m , 

2TOpj 

it can be derived easily that 

£ _ l-e(fc,a) 
* ~ l + e(k,a)'' 
fj.(k,a) = 1 + e(k, a). 



(194) 



l{k,a 



(195) 



These results are valid for all the models which can be 
described by a field tracking the minimum of the effective 
potential since before BBN. More precisely we find that 



fi(k,a) 

l{k,a) 



2^2 



{l + 2/3 2 )k 2 + m 2 a 

k 2 + m 2 a 2 
(l-2/3 2 )fc 2 +ro 2 a 2 
(l + 2/3 2 )fc 2 +m 2 a 2 



(196) 



These are closely related to the popular parameterisation 
of modified gravity used in the literature. Here they are 
valid for any model of modified gravity at the linear level 
of cosmological perturbations as long as the background 
cosmology is described by a scalar field slowly evolving in 
time and following the time dependent minimum of the 
effective potential where m 2 3> H 2 . 
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As a numerical illustration, in Fig. 5 we have com- 
pared the function j(a, k) calculated using three differ- 
ent methods: (1) the full numerical solution as shown by 
the black solid curve, (2) the value obtained by using the 
definitions in Eq. (192), the analytical approximation for 
5<j) in Eq. (184) and $jv solved from the Poisson equation 
numerically (the red dashed curve) and (3) Eq. (196) as 
shown by the blue dotted curve. We can see that the lat- 
ter two agree with each other very well, showing that the 
parameterisation given in Eq. (196) works very well in 
practice and describes the statistical properties of linear 
perturbations. 

The full numerical solution, however, again shows the 
oscillating behaviour, but the oscillation always centres 
around the averaged value defined by the previous for- 
mulae. As discussed earlier, over many oscillations there 
will be a cancellation and the net effect on a statistical 
observable today are the same for all the three curves. 



C. f(R) Gravity in the Jordan Frame 

Let us concentrate now on the case of f(R) gravity. 
The perturbations are then determined by 



H(k,a) 
l(k,a) 



|fc 2 + m 2 a 2 
k 2 + m 2 a 2 
|fc 2 + m 2 a 2 



4j,2 



k 2 + m 2 a 



2 n 1 



For large curvature models with m 
comes 



(197) 
moa~ r , this be- 



fi(k,a) 



j(k,a) 



4 k- a 3n+4 
3 raj. 



+ 1 



^ a 3«+4 + 1 ■■ 

+ 1 



2 k z „3n+4 

3 mi 



o rnf. 



(198) 



When n = |(r — 3) <C 1, we retrieve the phenomenolog- 
ical parameterisation [40] 



fi(k,a) 

-y(k,a) 






±]!LrA 



(199) 



3 1117. 



Our parameterisation Eq. (196) covers all the possible 
f(R) models. 



of the scalar field for scales k > 0.1 /iMpc -1 require full 
numerical simulations [53, 54]. 

The /3(a), m(a) parameterisation can completely spec- 
ify the nonlinear dynamics of </> with two temporal func- 
tions. Indeed, as we have seen above, one can reconstruct 
the potential V((f>) and the coupling function together 
with the background evolution <j)(a). Then one can study 
the non-linear evolution of the scalar field perturbation 
which, in the quasi-static limit, are governed by 



mpi mpi 



dV(<£) dV{4>) 



, (200) 



where the overbar means the background value. 

One can easily obtain dV((p)/d<f> analytically or numer- 
ically, and this can be used to solve the quasi-static dy- 
namics numerically. An advantage is that temporal func- 
tions 771(a), /3(a) completely specify the dynamics of (f>, in 
particular its spatial configuration, and there is no need 
for a /c-space paramctrisation. 

On linear scales, this is equivalent to the Jordan-frame 
description with the two spatially dependent function 
jj,(k, a) and j(k, a) being defined by e(k, a) which depends 
on the two functions 771(a) and /3(a), as given in Eq. (196). 
But in practice, working with two temporal functions is 
much more direct. Furthermore, the parameterisation de- 
scribed in Eq. (196) fails to faithfully describe the nonlin- 
ear effects or the environmental dependence. In essence, 
by going from 777(a) and /3(a) to fx(k, a) and j(k, a), one 
not only introduces spatial dependence but also loses the 
ability to describe nonlinear structure formation: in this 
sense, we may describe the approach using /j,(k, a) and 
"f(k, a) as the linear parameterisation of structure for- 
mation while 777(a) and /8(a) provide a fully non-linear 
parameterisation of modified gravity 11 . 

Past experience has shown that in modified gravity 
(e.g., chameleon and f(R)) models, nonlinear effects be- 
come important as soon as the linear perturbation result 
deviates from the corresponding ACDM prediction. This 
emphasises the importance of using full numerical simu- 
lations in the study of these models. However, the full nu- 
merical simulations arc generally very time and resource- 
consuming, and are therefore left for future work. 



VI. VARIATION OF CONSTANTS 

We have seen that the background evolution of the 
scalar field is specified by the time dependent mass and 
coupling functions. As the scalar field evolves, the particle 
masses and the gauge coupling constants change in time 



D. Non-linear Effects 

Matter clustering on galactic and cluster scales is an 
important probe of modified gravity. The nonlincarity in 
both the structure formation process and the dynamics 



Our parameterisation also provides a clear characterisation of 
the class of physical models (namely a scalar field coupled to 
matter) considered here, which is important in parametcrising 
modified gravity [48] , and not automatically incorporated in the 
(fi, 7) parameterisation. 
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too. The time variation of masses and gauge couplings 
is tightly constrained by laboratory experiments [58]. In 
this section, we analyse the time drift of the fine structure 
constant and the electron to proton mass ratio. 



A. The fine structure constant 

The scalar field also has an effect on gauge couplings 
and particle masses. The fermion masses are given by 



m F ((j)) = A((j>)m ha 



(201) 



where mbare is the bare mass in the Lagrangian. Mean- 
while, quantum effects such as the presence of heavy 
fermions lead to the potential coupling of </> to photons 

[59] 



S, 



1 



gauge 



4 5barc 



d /L xV=jB F (<j>)F llu F'" / , (202) 



where gbarc is the bare coupling constant and 



£f(0) = l + /3 7 — + 
mpi 



(203) 



The scalar coupling to the electromagnetic field would 
lead to a dependence of the fine structure constant on 4> 

as 



1 1 



a a ba 



■B F (</>), 



(204) 



The time evolution in the past is also particularly in- 
teresting. For symmetron models, we find that the time 
variation of a is 



Ha \mj aVl-tv) 3 



(207) 



Here, the time variation of a increases as one reaches 
the transition a*. This is a large variation which may 
happen in the recent past of the Universe and may have 
observable consequences in the emission lines of distant 
objects. 

It should however be noted that even though a/a can 
be very large, the relative difference of a between the 
earth and some other sparser place in the Universe is 
constrained to be less than 



Aa 



< 



TOP] 



3/3*/3 7 fi m0 ( — ) — (208) 

"I* / PmO 



If we instead consider a quadratic coupling to photons, 



B F ((t)) = 1 + 4f-0 2 , we find 
Aa 



< Al^- = 3A^ T fi m0 ( — 



Hp 
m 



P* 

PmO 



(209) 



where /3 7 = fiA\/ A2 = cfi^mpiA^. 

Interestingly, for both cases and for our fiducial param- 
eter values m* ~ 10 3 i?o, p* ~ p m o and /3 ~ /3 7 = 0(1) 
this term is of the same order as the claimed variation of 
a reported in [61]. 



implying that 



Q 

a 



-& 



mpi 



(205) 



where we have assumed that /3 7 ^>/mpi <C 1. Using the 
evolution equation we find that 



-£- « -9/3 7 /?fi m — 
Ha m 



(206) 



Hence the negative variation of the fine structure con- 
stant in one Hubble time is related to the small ra- 
tio H/m -C 1 and the couplings of <f> to matter and 
photons. The best experimental bound on the variation 
of a now comes from Aluminium and Mercury single- 
ion clocks [60]: || = (-1.6 ± 2.3) • lO^V -1 - Tak- 



ing H 



1.5 ■ 10 yr, we get the conservative bound 



l;§M < 2 • 10 7 . As a result, the experimental bounds 



on the time variation of a lead to constraints on /3o/3 7 o 

as /3 /3 7 o < 0.8 • lO" 7 ^!. For models with f3 = (9(1), 

f^roO ~ 0.25 and mo/Ho « 10 3 where effects on large 
scale structure are present, /3 7 o < 0.1, which is a much 
tighter bound than present experimental ones /3 7 o < 10 11 

[55]. 



B. The variation of masses 

Fundamental fermions such as the electrons have a uni- 
versal mass dependence m F = A(^)mbare> implying that 



m F ~ Q 2c- H 



Hm F 



9^Q r 



,2- 



(210) 



Nucleons such as the proton have a mass given by the 
phcnomenological formula 

m p = CqcdAqcd + b u m u + b d m d + C p a, (211) 



where Aq C d ~ 217MeV is the QCD scale, b u + b d ~ 6, 
b u - b d ~ 0.5, Cqcd ~ 5.2, mJr° ~ 5MeV, m h d aie ~ 
lOMeV and C p a\, aro ~ 0.62MeV. Assuming conserva- 
tively that Aqcd is scalar independent, we get 



Hm 



QQm pz_ ( »»">* 2 Ud " bd a - ^^/3 7 



H 2 _ ( b a ml"~ + o d m d - 



m 



111, 
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It is particularly important to study the variation of 

nip 



111 



(212) 
(213) 



p 
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from which we find that its time variation is positive for 
modified gravity models: 



9Q m/ 0- 



M 2 



p+ C p a haic ^ 



(214) 



The current experimental constraint is -|o = (—3.8 ± 
5.6)10 _14 yr _1 which yields the upper bound on /3o: /3 2 < 
10" 5 ^|. For /3 = Oil), this entails that m Q /H a > 10 25 . 
Again for symmctron models, the electron to proton mass 
ratio would vary rapidly in time around the transition 
time a*. It would be interesting to study if such a varia- 
tion could have relevant effects on the physics of distant 
objects. 



VII. CONCLUSION 

We have developed a novel parametrisation of modified 
gravity models first presented in [5]. Starting with the 
time-evolution of the mass and the matter coupling of a 
scalar field in the cosmological background, we have been 
able to reverse engineer the complete dynamics of these 
models in a simple way. 

We have applied these results to well-known modified 
gravity models: chameleons, f(R) gravity, dilatons and 
symmetrons. In each case, we have explicitly given the 
mapping and the full reconstruction. We have also shown 
how one can apply local constraints using this formalism 
and then use it to make predictions for linear cosmolog- 
ical perturbations. 

New classes of models can be engineered in a more 
intuitive way than starting from a Lagrangian. The La- 
grangian itself can be completely reconstructed. One only 
needs to specify two functions whose physical meaning is 
easily grasped: namely the mass (the inverse range of the 
fifth-force) and the coupling to matter. 

The real strength of this approach compared to ex- 
isting parameterisations in the literature is that we can 
reconstruct the whole theory at the linear and non-linear 
levels and be sure that it corresponds to a concrete 
physical model defined via a Lagrangian. This effectively 
superscd existing parameterisations of modified gravity 
with a screening mechanism by being able to make pre- 
dictions for non-linear clustering of matter via iV-body 
simulations. This will be the subject of future work. 
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Appendix A: The Cavity Constraint 

In this appendix, we will explicitly develop the calcu- 
lation leading to the cavity constraint for chameleon and 
f(R) models. 

Consider a cavity of radius L with a residual density 
Pcav *C p c where p c is the density of the bore surrounding 
the cavity. The field inside the cavity is <^ cav and deviates 
slightly from this value across the cavity. Expanding the 
effective potential around 4> cav and putting Scf> = 4> — cav , 
we have 

1 d 
r 2 dr 



-o X~ r T 5 ^ ~ m cav^ = V cS,4> (0cav) ( A1 ) 



where m cav is the scalar field mass inside the cavity and 
V^ff.0 = dV e fi(<f>)/d(f) is nonzero unless cav minimises the 
effective potential. Inside the cavity the solution is 



= A 



sinh(m cav r) V cS ^ (</> cav ) 



outside the cavity we have 



B- 



(A2) 



(A3) 



where A, B are constants of integral, <j) c is the minimum 
of the effective potential outside the cavity and m c the 
mass at that minimum. Matching at r = L, we find that 



a m c L 



B 



1 + m c L 



and 



A 



1 + m c L 



[sinh (m cav L) - m cav L] A, 



sinh (m cav L) 



(A4) 



1 + m r L 



V cS , 



(f> \(Pca,v 



in- 



(A5) 



Evaluating the solution at the origin and putting 6<f>(r = 
0) = we have 



A 



V c 



efi,<f> Wcav 



111' 



This leads to 



1 



sinh {m cav L) 



4a V m 



2 

cav 



Kff,0 (^cav) 



(A6) 



(A7) 



where we have used m c L ^$> 1. 

For potentials V ~ 1/0" and as long as 4> cav is much 
less than the effective minimum in the cavity we have 

sinh (m cav £) 

r = n : ( A8 ) 

which implies that 

m cav L = (9(1), (A9) 

where m cav is dominated by the potential term. 
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